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ABSTRACT 

Certain  ORION  missions  may  require  three  axis  stabilization.  Since  ORION's 
physical  size  severely  limits  its  onboard  fuel  storage  capability,  passive  stabilization 
techniques  warrant  investigation.  This  paper  shows  the  development  of  linearized 
equations  of  motion  and  regions  of  stability  with  respect  to  gravity  gradient  stabiliza- 
tion. Gravity  gradient  stabilization  by  itself  provides  little  yaw  restoring  torque:  there- 
fore, additional  torque  generating  devices  are  necessary  to  augment  the  gravity  gradient 
effect.  Control  moment  gyros,  reaction  wheels,  and  magnetic  torquers  will  be  investi- 
gated as  to  their  suitabilitv  for  ORION. 
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I.     INTRODUCTION 

A.     BACKGROUND  ON  ORION 

The  philosophy  behind  the  design  and  development  of  the  Orion  satellite  is  being 
driven  in  part  by  economic  considerations.  Small,  affordable,  general  purpose  satellites 
must  be  developed  without  sacrificing  reliability.  The  preliminary  groundwork  on  the 
feasibility  of  the  ORION  concept  and  design  is  detailed  in  [Ref.  1]. 

Orion  is  currently  being  designed  to  be  launched  from  a  Get  Away  Special  (GAS) 
cannister  aboard  the  Space  Shuttle.  The  GAS  cannister  physically  limits  the  size  of 
ORION  to  approximately  35"  in  height  and  19"  in  diameter.  Figure  1  shows  the  Get- 
Away  Special  Cannister  attached  to  the  side  of  the  shuttle  bay. 
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Figure  1.  NASA  Get-A>vay-Special  (GAS)  Cannister  [Ref.  1,  p.  30] 


The  Get-Away-Special  (GAS)  Cannister  is  designed  to  launch  its  satellite 
autonomously  on  command.  A  cross  sectional  view  of  the  (GAS)  Cannister  is  shown  in 
Figure  2  on  page  2. 

ORION  should  be  able  to  meet  the  requirements  of  standard  or  typical  payloads. 
Some  of  these  requirements  are  listed  by  Reference  2  and  given  in  Figure  3  on  page  3. 
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Figure  2.       Cross  section  of  GAS  Cannister  [Ref.  1,  p.  42] 


1.     Background  on  attitude  control 

The  success  of  any  satellite  is  dependent  on  its  attitude  control  system.  To  be 
mission  capable,  the  vehicle  must  first  acquire  a  desired  or  specified  orientation  in  space, 
and  then  maintain  this  orientaion  within  mission  limits  given  such  external  factors  as 
magnetic  fields,  aerodynamic  drag,  solar  pressure,  gravitational  effects,  and  other  dis- 
turbances. Deviations  from  desired  orientations  are  detected  through  the  use  of  sensors. 
Considerations  for  an  attitude  control  system  include: 


(1)  pointing  accuracy 

(2)  fuel  consumption 

(3)  thermal 

(4)  power  constraints 

(5)  satellite  lifetime 

(6)  response  time 


Mass 

3  2  lbm 

Volume 

2.3  6  ft3 

Power  ' 

34  watts 

DataRate 

5000  bits/sec 

Orbit 

200-000  nm  circular 

Inclination 

0'-30'  or  60e-120° 

Instruments 

Particle  counter  or 

Lens 

Figure  3.      Summary  of  typical  Navy/STP  pay  load  requirements 


The  weight  given  to  each  of  these  considerations  ultimately  lies  in  the  specific  mission 
requirements  designated  for  each  satellite.  The  direction  of  this  thesis  will  be  to  inves- 
tigate the  feasibility  of  achieving  sufficient  performance  utilizing  a  gravity  gradient  con- 
trol system.  The  standard  of  +/-  1  degree  for  each  axis  will  be  used.  This  criteria  is 
sufficient  to  enable  ORION  to  satisfy  most  SIP  (space  test  program)  mission 
requirements.! Ref.  1,  p.  93) 

1  o  achieve  stabilization,  many  methods  of  generating  controlling  torques  have 
been  developed.  These  methods  are  divided  into  either  active  or  passive  categories.  Each 
of  these  categories  has  their  relative  advantages  and  disadvantages  [Figure  4  on  page  4 
and  Figure  5  on  page  4]. 

In  this  study,  the  primary  focus  is  on  a  gravity  gradient  stabilization  augmented 
by  other  torque  generating  devices  to  achieve  the  desired  +  /-  1  degree  goal  for  each  axis. 
A  discussion  of  the  equations  of  motion  and  gravity  gradient  stabilization  theory  is 
presented  in  Chapter  2.  Disturbance  torques  and  other  torque  generating  devices  are 
discussed  in  Chapter  3.  Chapter  4  presents  the  results  of  stabilization  schemes  derived 
in  Chapters  2  and  3.  Observations  and  conclusions  are  discussed  in  Chapter  5. 
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Figure  4.       Passive  Stabilization  Methods  [Ref.  4,  p.  19] 
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Active  Stabilization  Methods  [Ref.  4,  p.  19] 


II.     GRAVITY  GRADIENT  STABILIZATION  THEORY 


A.     COORDINATE  DEFINITION 

Before  beginning  any  discussion  on  attitude  control  a  means  of  defining  the  coordi- 
nate system  is  necessary.  The  standardized  coordinate  axes  for  gravity  gradient  systems 
are  as  follows:  the  x  axis  points  in  the  direction  of  satellite  motion,  the  y  axis  is  normal 
to  the  orbit  plane,  and  the  z  axis  points  towards  earth's  geocenter.  These  axes  corre- 
spond to  the  roll,  pitch  and  yaw  axes,  respectively.  Figure  6  from  Reference  3  defines 
roll  (</)),  pitch  (#),  and  yaw  {\[/)  attitude  errors  as  well  as  showing  the  standardized  co- 
ordinate svstem. 


COORDINATE  DEFINITION 
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Figure  6.      Coordinate  Definition  for  Gravity  Gradient  Systems 


B.     COORDINATE  TRANSFORMATION 

The  relationship  between  two  coordinate  systems  is  essential  for  attitude  control. 
There  are  methods  available  to  describe  this  relationship  of  which  directional  cosines, 
quaternions,  and  Euler  angles  will  be  briefly  discussed.  Figure  7  on  page  6  shows  various 
methods  of  coordinate  transformations  along  with  their  respective  advantages  and  dis- 
advantages. 


PARAMETERIZATION 

NOTATION 

AOVANTACES 

DISADVANTAGES 

COMMON  APPLICAllONS 

DIRECTION 
COSINE 
MATRIX 

A  -  |A    | 

NO  SINGULARITIES 
NO  TRIGONOMETHIC  FUNCTIONS 
CONVENIENT  PRODUCT  RULE 
FOR  SUCCESSIVE  ROTATIONS 

SIX  REDUNOANT  PARAMETERS 

IN  ANALYSIS    TO  TRANSFORM 
VECTOHS  FROM  ONE  REFER. 
ENCE  FRAME  TO  ANOTHER 

tULER 

AXIS/ANCLE 

«.  » 

CLEAR  PHYSICAL  INTERPRETATION 

ONE  REDUNDANT  PARAMETER 
AXIS  UNDEFINEO  V»HEN  SIN  ••0 
TRIGONOMETRIC  FUNCTIONS 

COMMANDING  SLEW  MANEUVERS 

EULER 

symmetric 
Parameters 

iouaternioni 

Id 

NO  SINGULARITIES 
NO  TRIGONOMETRIC  FUNCTIONS 
CONVENIENT  PRODUCT  RULE   FOR 
SUCCESSIVE  ROTATIONS 

ONE  REOUNDANT  PARAMETER 
NO  OBVIOUS  PHYSICAL  INT  ER- 
PRETATION 

ONBOARD  INERTIAL  navigation 

cises 

VECTOR 

I 

NO  REDUNDANT  PARAMETERS 
NO  TRIGONOMETRIC  FUNCTIONS 
CONVENIENT  PRODUCT  RULE  FOR 
SUCCESSIVE  ROTATIONS 

INFINITE  FOR    IBP  DEC  ROTA- 
TION 

ANALYTIC  STUOIES 

EULER 

ANCLES 

».  ».  * 

HO  REDUNDANT  PARAMETERS 
PHYSICAL    INTERPRET  ATION  IS 
CLEAR  IN  SOME  CASES 

TRIGONOMETRIC  FUNCTIONS 
SINGULARITY   AT  SOME   » 
NO  CONVENIENT  PRODUCT  RULE 
FOR  SUCCESSIVE  ROTATIONS 

ANALYTIC  STUDIES 
INPUT /0U1PUT 
1  ONBOARD  ATTITUDE  CONTROL  OF 
1  AXIS  STABILIZED  SPACECRAFT 

Figure  7. 


Coordinate  Transformation  Methods  [Ref.  4,  p.  412) 


1.     Direction  cosines  method 

The  directional  cosine  method  involves  a  3  x  3  rotation  matrix  A„  where 


"^11 

'4\2     ^13 

4/  = 

^21 

/J22     ^23 

= 

A3} 

AU     ^33_ 

X*x  X»y  X*z 
Y'X  Y*y  Y.z 
Z  •  x    Z  *y   Z  •  z 


(1) 


where 


X  •  x  =  cos  6  cos  iff 

X»y  =  cos  \\i  sin  0  sin  4>  +  sin  ^  cos  </> 

X » z  —  sin  ip  sin  <fi  —  cos  \}/  sin  6  cos  4> 

Y  *x  =  —  sin  \p  cos  6 

Y  »y  =  cos  t/>  cos  if/  —  sin  \}/  sin  6  sin  <p 
Y » z  =  sin  \p  sin  6  cos  <f)  +  cos  \[/  sin  0 
Z  •  x  =  sin  6 

Z  »y=  —  cos  6  sin  4> 
Z  •  z  =  cos  0  cos  6 


2.     Quaternions 

The  quaternion  method  or  as  it  is  sometimes  called,  the  Euler  Symmetric  Pa- 
rameter method,  uses  the  Euler  axis  and  rotation  angle  (a  single  rotation  angle  a  about 
a  single  axis  E)  method  to  describe  the  relationship  between  two  coordinate  systems 
[Figure  8  on  page  7J. 


Figure  8.      Rotation  about  the  Euler  Axis 

The  quaternion  is  a  four  component  vector  containing  the  same  information  as  the  Euler 
angle/ axis  transformation  in  which  one  coordinate  axis  is  related  to  another  and  is  of  the 
form 


<7  =  Ui'»  <7iA  ft*>  ?J 


(3) 


where 


<7,  =  £,siny 
q2  =  Ey  sin  — 

r  « 

ft  =  Ez  S1«  y 

a 
<74  =  cos  y 


The  symbols  i,  j,  and  k  satisfy  the  following  conditions  [Ref.  4,  p.  758J 


.2         .2         ,2         , 

i   =J   =  k   =  1 
ij  =  -ji  =  k 

jk  =  -kj=  i 
ki  =  -ik  =j 


and 


2  2  2  2 


1=1 


When  used  for  control  purposes,  [Ref  5]  defines  the  first  three  elements  of  the 
quaternion  {qu  q2,  g3)  as  the  respective  roll,  pitch,  and  yaw  rotational  errors: 


LErotl>  epttch>  €ya*~\  =  ^C^i,  ^2-  ^3^ 

=  2   Ex  sin  -f-  ,  Ey  sin  —  ,  Ez  sin  — 
?z\_y.Ex,  o>.Ev,  o.E2~\ 


(7) 


The  primary  advantage  of  the  quaternion  method  is  that  the  computational  time 
can  be  reduced  by  more  that  40%  over  the  equivalent  directional  cosine  matrix  solution 
in  operations  that  require  successive  coordinate  system  rotations  [Ref  5].  The  transfor- 
mation from  coordinate  system  A  to  coordinate  system  B  can  be  described  by  the 
quaternion  q  =  \_qxi,  qj,  q3k,  q^.  and  the  second  transformation  from  system  B  to  C  can 
be  described  by  q'  =  \_q\i.  q'J.  q\k,  q\~\.  The  transformation  from  system  A  to  system 
C  may  be  described  by  the  quaternion  q"  =  \_q" \i,  q"J,  q"3,  q'\~]  where 


q"  =  qq'  =  (qxi  +  qj  +  q2k  +  qM  xi  +  q'J  +  q'3k  +  q'4) 


(8) 


Multiplying  out  and  substituting  in  the  conditions  stated  above,  the  matrix  form  be- 
comes 


q  4  q  3  -?2  fi 

-q's  q\  q\  q\ 

q'i  -q'\  q\  4* 

-<f\  -q\  -q'i  q\ 


ft 

<?3 
?4 


(9) 


Because  of  the  quaternion  property  of  interchangeability,  a  total  of  16  multiplications 
and  12  additions  are  required  to  completely  define  the  quaternion  describing  the  trans- 
formation from  coordinate  system  A  to  coordinate  system  C.  To  accomplish  the  same 


transformation  using  directional  cosines  would  require  multiplying  two  3x3  matrices, 
involving  27  multiplications  and  18  additions.  The  computational  time  savings  for  suc- 
cessive rotations  transformations  becomes  apparent.  The  major  drawback  for  the  use 
of  quaternions  is  that  the  numbers  of  the  quaternion  do  not  physically  represent  the 
transformation  from  one  system  to  another. 
3.     Euler  angles  method 

The  Euler  angles  utilize  three  different  rotations  (angles),  defining  the  orien- 
tation of  a  body  with  respect  to  an  inertial  reference  frame.  These  angles  are  defined  as 
$  (about  the  x  axis),  6  (about  the  y  axis),  and  </>  (about  the  z  axis).  These  axes,  by  es- 
tablished convention,  may  also  be  referred  to  as  1  (for  x),  2  (for  y),  and  3  (for  z)  axes, 
respectively. 

The  321  Euler  transformation  described  below  starts  with  the  inertial  coordinate 
system  XYZ.  Figure  9  on  page  10. a.  shows  the  first  rotation  about  the  Z  axis  through 
an  angle  \p,  which  produces  xy\Z.  Figure  9  on  page  lO.b.  shows  the  rotation  of  6  about 
y\  ,  producing  xyxzv  Figure  9  on  page  10. c,  showing  the  rotation  of  <j>  about  x,  produces 
xyz.  The  rotation  described  in  Figure  9  on  page  10. a.  about  the  z  axis  produces 
Figure  10  on  page  10.   The  transformation  from  XYZ  to  xy\Z  is 

jCj  =  X  cos  \\j  +  }'  sin  i/> 
yx  =  —X  sin  (//  +  Y  cos  \]/ 
Z  =  Z 

or  written  in  matrix  form 


z 


cos  \jj      sin  \j/    0 

—  sin  \p    cos  \p    0 

0  0        1 


X 
Y 
Z 


(11) 


figure  9.      Euler  Angles 


Figure   10.      Rotation  about  the  Z  axis  through  the  Angle  t/> 


Writing  the  transformations  of  Figures  9.b.  and  9.c.  in  matrix  form  yields 


X 

cos0   0 

—  sin  6 

x\ 

y\ 

= 

0       1 

0 

y\ 

Z| 

sin0    0 

cos  0 

z 

(12) 
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and 


X 

1 

0 

0 

X 

y 

= 

0 

COS  0 

sin  0 

y\ 

z  _ 

0 

—  sin  0 

cos  0 

ji 

(13) 


Combining  eqs(l  1-13)  into  one  equation  yields  the  complete  321  coordinate  transfor- 
mation from  the  inertial  reference  frame  to  the  body  reference  frame 


1         0  0 

0     cos  0      sin  0 
0    —  sin  0    cos  0 


cos  6   0    —  sin  $ 

0       1         0 
sin  6    0     cos  6 


cos  \p      sin  \p    0 

—  sin  \j/    cos  \j/    0 

0  0       1 


/ 
J 

K 


(14) 

cos  6  cos  41  cos  $  sin  ^  —  sin  6 

-  cos  0  sin  0  +  sin  0  sin  0  cos  ij/      cos  0  cos  \\)  +  sin  0  sin  6  sin  t//      sin  0  cos  6 
sin  </>  sin  \j/  +  cos  0  sin  6  cos  <//      —  sin  0  cos  \jj  +  cos  0  sin  6  sin  ^    cos  0  cos  6 


This  is  just  one  of  the  12  different  Euler  coordinate  transformations  schemes  available. 
Figure  11  on  page  12  lists  all  possible  Euler  angle  combinations  and  their  resulting 
transformation  matrices.  It  should  be  noted  that  eq(14)  does  not  agree  with  the  321 
Euler  transformation  matrix  of  Figure  11  on  page  12.  The  321  matrix  of  Figure  11  on 
page  12  contains  a  error  that  interchanges  <f>  with  \j/  and  \J/  with  0. 


X 
Y 
Z 
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i  >  r  r  .  i 

rin  r m  n<ci  t 

nuntstNiAnon 


I  •  f  -  3 


1    J    I 


J    l    J 


J    11 


MAIniX  A 


»vrt - 1 

rin  tn  />nr.i  s 

ntpnEsEMiAiion 


i    j    I 


i- J    I 


i   i   ; 


j  i  j 
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-stcr 
w 

c*ci 

-SI 
S*Cl 

c»c« 
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II 
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-Sf 
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-S*C*'  C*S*S* 
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c« 

S"S-» 
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Figure  11.      Table  of  Euler  angle  combinations  [Ref.  4,  p.  764) 
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C.     RIGID  BODY  EQUATIONS 

The  general  motion  of  a  particle  and  its  corresponding  momentum  have  been  written 
about  extensively  in  various  publications.  A  summary  of  the  relevent  derivations  and 
equations  of  motion  from  [Ref  6]  are  presented  here. 

1.     Derivative  of  a  Vector  in  a  Rotating  Coordinate  System 

The  time  derivative  of  a  vector  in  a  rotating  coordinate  system  consists  of  the 
rate  of  change  of  the  vector  relative  to  the  moving  axis  and  the  rate  of  change  due  to  the 
rotation  of  the  axis.  Figure  12  shows  a  vector  in  a  rotating  system  where  XYZ  represent 
an  inertial  system  and  xyz  are  the  set  of  axes  rotating  with  angular  velocity  co  relative 
to  XYZ.  If  unit  vectors  along  the  X,  Y,  and  Z  axes  are  i,j,  and  k,  respectively,  the  vector 
r  can  be  written  in  the  form 


r  =  xi  +  yj  +  zk 


(15) 


Figure  12. 


Vector  in  a  rotating  system  (Ref.  6,  p.  107] 


The  time  derivative  of  r  is 


r  —  xi  4-  xi  +  yj  +  yj  +  zk  +  zk 


(16) 


From  [Ref.  7],  it  has  been  shown  that  the  time  derivatives  of  the  unit  vectors  are 


/  =  co  x  / 


(17) 
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J  =wxj 


(18) 


k  =  wxk  (19) 

Defining  rtel  as  the  time  rate  of  change  of  r  relative  to  xyz,  eq(16)  becomes 

r  =  rrel  +  co  x  r  (20) 

2.     General  Particle  Motion 

Translational  motion  of  a  particle  with  references  to  two  coordinate  systems  is 
shown  in  Figure  13. 
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Figure   13.      Translational  motion  of  a  particle  [Ref.  6,  p.  107] 


The  position  vector  of  a  particle  relative  to  coordinate  system  XYZ  is 


R  =  Rc  +  r 


(21) 


where  #c  is  the  vector  from  the  axis  of  the  XYZ  coordinate  frame  to  the  xyz  axis  frame 
and  r  is  the  differential  (position)  vector  in  the  xyz  frame.  The  velocity  of  the  particle 
is  the  derivative  of  eq(21)  with  respect  to  time 


V=R  =  Rc  +  r 


(22) 
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The  velocity  at  the  origin  Rc  =  Vc  combined  with  eq(20)  yields 

V=  Vc+  Vrel+coxr  (23) 

3.     Momentum  of  a  Rigid  Body 

A  particle's  momentum  is  the  product  of  its  velocity  and  mass.    The  linear  mo- 
mentum of  the  particle  of  Figure  13  on  page  14  is 


=  j  Vdm  (24) 


P 

'm 


where  V  is  the  particle's  velocity  and  m  is  its  mass.  Recalling  that  Vrel  is  the  time  rate 
of  change  of  r  relative  to  the  xyz  coordinate  system,  and  assuming  that  the  xyz  axis  is 
centered  in  the  rigid  body,  Vrel  =  0.    Substituting  eq(23)  into  eq(24)  with  V„t=0, 

P  =      (Vc  +  to  x  r)  dm  =  mVc  +  on  x  \   r  dm  (25) 

By  choosing  the  center  of  mass  of  the  rigid  body  as  the  center  of  the  rotating  coordinate 
system,  f  r  dm  =  0  by  definition.  This  simplifies  eq(25)  to 

P  =  mVc  (26) 

where  P  is  the  linear  momentum  of  the  rigid  body,  m  is  its  mass,  and  Vc  is  the  velocity 
of  its  center  of  mass. 

The  angular  momentum  of  a  particle  dm  about  point  C  is  defined  as  the  moment 
about  C  of  the  linear  momentum  given  by 

hc  =  rxVdm  (27) 

Substituting  for  V  yields 


Hc  =      r  x  V  dm  =      r  x  ( Vc  +  co  x  r)  dm 


(28) 


Rearranging  terms 
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Hc  = 


r dm)  x  Vc  + 


r  x  (co  x  r)  dm 


(29) 


By  choosing  the  center  of  mass  of  the  rigid  body  as  the  origin  of  the  rotating  coordinate 
system,  eq(29)  reduces  to 


Hc  =      r  x  (co  x  r)  dm 


(30) 


Reducing  the  angular  velocity  vector  o>  into  its  components 


to  =  coxi  +  cOyj  +  cozk 


(3D 


Combining  equations(15),  (30),  and  (31)  yields 


+    [    -IX;C0X    ~    Iy2(0y    +  /„«J* 


(32) 


where 


'xx 


(y   +  z  )  </m 


(33) 


Iyy  =  \    (x2  +  z2)  dm 


(34) 


':; 


dm 


(35) 


are  the  moments  of  inertia  about  the  x,y,  and  z  axes,  respectively,  and 


-i 


Jjcv  =       *>'  dm 


(36) 
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-L- 


Ixz  =  I   .rz  dm 


(37) 


l 


/y2  =  I  yz  dm 


(38) 


are  the  products  of  inertia.  Equation(32)  written  in  matrix  form  becomes 


Hy 

H2 


I„     -/„    -/ 


xx  *xy  *xz 


*xy       *yy 

-L,    -l 


lyz 


xz         'yz         *zz 


W, 


CO, 


La)2 


(39) 


4.      Equations  of  motion/ Euler's  Moment  Equations 

The  derivation  of  the  following  Euler's  moment  equations  are  summarized  from 
[Ref  6].   Recalling  that  angular  momentum  hc  relative  to  the  point  C  is 


hr=  r  x  m\ 


(40) 


its  derivative  with  respect  to  time  is 


hc=  r  x  mv  +  r  x  mr 


(41) 


The  moment  of  the  force  about  a  point  C  by  definition  is 


Mc  =  r  x  F 


(42) 


Recall  Newton's  second  law  describing  force: 


F '  =  ma  =  mv 


(43) 


Substituting  eqs(22)  and  (43)  into  eq(41)  yields 


hc  =  r  x  m(Rc  +  r)  +  r  x  F 


(44) 


Substituting  eq(42)  into  eq(44)  and  rearranging  terms  yields  the  moment  equation 
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MC  =  HC  + 


(Rc  x  f)  dm 


(45) 


where  Hc  is  the  angular  momentum  of  the  rigid  body  about  point  C.  From  eq(20)  the 
latter  part  of  eq(45)  can  be  written  as 


Rcx  r  dm  =  Rcx  I  rrel  + 


{axirdm)) 


(46) 


Noting  that  point  C  is  the  center  of  mass  f  r  dm  =  0  and  rrtl  =  0  for  a  rigid  body,  eq(46) 
reduces  to 


Mc  =  Hc 


(47) 


The  moment  of  external  forces  about  the  center  of  mass  of  a  rigid  body  is  equal  to  the 
time  rate  of  change  of  the  angular  momentum  of  the  body  about  the  center  of  mass  [Ref. 
7].    In  terms  of  a  rotating  axis,  the  moment  equation  of  eq(47)  becomes 


M,  =  H,rpl  +  cd  x  Hr 


lc       ilcrel 


(48) 


Substituting  eq(39)  in  for  H,  the  components  of  the  moment  equation  can  be  written  as 


Mx  =  Hx  +  ooyHz  —  (ozHy 


(49) 


My  =  Hy  +  oozHx  —  (oxHz 


(50) 


Mz  =  Hz  +  (oxHy  —  coyH: 


y"x 


(51) 


where  Mx^  are  the  moment  equations.  In  this  study,  an  assumption  is  made  that  the 
principle  axes  of  the  moment  of  inertia  are  the  rotating  axes  x,y,z;  therefore,  the  pro- 
ducts of  inertia  (Iiy,  /„,  Iy2)  are  equal  to  zero.  The  angular  momentum  components  re- 
duce to 


Hx  =  Ixx^x 


(52) 
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Hy    =    IyyCOy 


(53) 


Hz  -  I22co2 


(54) 


Substituting  eqs(52-54)  into  the  moment  equations  yields  Euler's  moment  equations 

MX     =      lXX(bX     +     0Jy0J2{lZ2     -      iyy)  (55) 


My    =     IyyUy    +     fOjfiijJ^    ~     I„) 


M2     =      l22l02     +      0}X(Oy{Iyy     ~      /jj 


(56) 


(57) 


Referring  to  the  coordinate  system  shown  in  Figure  6  on  page  5,  the  angular 
velocities  in  terms  of  orbital  rate  co0  and  attitude  error  angles  pitch(0),  roll(c/>),  and  yaw 
(\p)  from  [Ref.  6,  p.  130]  are 


(O, 


CO, 


LW2J 


1         0  -  sin  0 

0     cos  c/>      cos  6?  sin  <j) 
0    —  sin  c6    cos  0  cos  c/> 


-   Or 


cos  0  sin  \j/ 

cos  0  cos  t//  +  sin  4>  sin  0  sin  \j/ 
—  sin  c/>  cos  \p  +  cos  c/>  sin  0  sin  ^ 


(58) 


When  using  small  angle  approximations,  eq(58)  reduces  to 


cox 

4>  -  (Oq\I/ 

(Oy 

= 

e-co0 

_WZ_ 

<A  +  ^0^ 

(59) 


Note  that  the  orbital  rate  is  defined  as 


COn  = 


R 


(60) 


where  ix  is  the  gravitational  constant  and  R  is  the  radius  of  the  earth  plus  the  satellite's 
altitude.  For  the  purpose  of  station  keeping,  the  small  angle  approximations  are  valid. 
Substituting  eq(59)  into  the  rigid  body  moment  eqs(55-57)  yields  the  linearized  rigid 
body  dynamic  equations 
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MxRB  =  (xx0  +  wo(A>v  ~  hz)4>  ~  W0(7xx  +  hz  ~  hy)$ 


yy ' 


(61) 


MyRB  —  lyyQ 


(62) 


M2RB  =  hz^  +  W0(/Vv  -  ^xx)"A  ~  Wo(/„  +  lZZ  ~  Iyy)4> 


yy) 


(63) 


D.     GRAVITY  GRADIENT  EFFECT 

Attitude  control  using  gravity  gradient  stabilization  works  on  the  principle  that 
given  an  object  with  an  asymmetrically  distributed  mass,  the  object  will  tend  to  align  it- 
self with  its  minimum  moment  ofinertia  axis  along  the  local  vertical  towords  the  earth. 

The  force  due  to  gravity  on  a  mass  element  is 


dF  = 


y.edm 


R 


(64) 


where 


R  =  Rt>  +  r 


(65) 


and 


f*t 


GMe  =  geR; 


(66) 


where  G  =  universal  gravity  constant,  Me=  the  mass  of  the  earth,  ge=  gravity  acceler- 
ation at  the  earth's  surface,  and  Re  =  the  radius  of  the  earth.  Reference  8  shows  the 
torque  generated  about  the  mass  center  of  the  object  is 


dM=rxdF=- 


Hedm 
\xedm 


rx  R 


rx(RQ  +  r) 


R 

Hgdm 

3     rx  Rq 
R3 


(67) 


Cubing  eq(65)  yields 


K3  =  *„3 


3rfip        3r2^         r3 

9  "^  1 


Rn 


0 


R 


0 


R 


0 


(68) 
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Noting  that  r4R0,  tne  higher  order  terms  drop  out.   Inversion  of  eq(68)  yields 

j  J     /        3r  •  Rr 


/T         Rr 


1 


R£ 


(69) 


Substituting  eq(69)  into  eq(67)  and  integrating  to  find  the  moment  yields 


M. 


=  -He\   rdmxR0  +  —f- 

Jm  An     Jr 


3ne   f    (rxfloXr.Ro) 


rdm  x  £0  +  -^f ^ —  rfm  (70) 


Given  that  the  center  of  mass  is  the  origin  of  the  body  frame,  f  rdm  =  0,  eq(70)  reduces 
to 


M« 


/<n    Jm 


Xr  •  «o)  dm 


(71) 


In  the  reference  system  defined  by  Figure  14,  the  K  (Z)  vector  is  pointing  towards 
the  earth;  therefore,  the  vector  R0  in  the  reference  coordinate  frame  is 


R0r=RoK 


(72) 


Figure   14.      Gravity  gradient  torque  [Ref.  6,  p.  1131 


21 


Restating  the  Euler  321  rotation  matrix  described  earlier 


1         0  0 

0     cos</>      sin  0 
0    —  sin  4>    cos  </> 


cos  0    0    -  sin  0 

0       1         0 
sin  0    0     cos  0 


cos  i>      sin  i//    0 

—  sin  \jf    cos  i//    0 

0  0        1 


/ 
J 
K 


(73) 

cos  0  cos  \jj  cos  0  sin  ^  —  sin  0 

-  cos  4>  sin  (/)  +  sin  </>  sin  0  cos  ^      cos  </>  cos  i/>  +  sin  <f>  sin  0  sin  \\i      sin  </>  cos  0 
sin  </>  sin  t/'  +  cos  0  sin  0  cos  t//      —  sin  </>  cos  t/'  +  cos  4>  sin  0  sin  i/^    cos  </>  cos  0 


7 
J 
K 


Transforming  R0r  into  the  spacecraft's  fixed  coordinates  i.  j,  and  k: 


&0r  = 


cos  0  cos  \j/  cos  0  sin  </f  —  sin  0 

-  cos  (j)  sin  0  +  sin  </>  sin  0  cos  i//      cos  </>  cos  </>  +  sin  <£  sin  0  sin  \p      sin  0  cos  0 
sin  4>  sin  i//  +  cos  0  sin  0  cos  \\i      —  sin  (/>  cos  \jj  +  cos  </>  sin  0  sin  i//    cos  </>  cos  0 


0 
0 
A' 


which  reduces  to 


R0r  =  R0{  —  sin  0/  +  sin  0  cos  0/  +  cos  </>  cos  0A-) 


(75) 


Substituting  eq(75)  into  eq(71)  yields  the  gravity  gradient  moment  equations 


A/„  = 


ft? 


(/,,  —  Iyy)  sin  0  cos  0  cos^O 

—  (Ixx  —  Iz2)  sin  0  cos  0  cos  </> 

—  (/jy,  —  7^)  sin  0  cos  0  sin  4> 


(76) 


Replacing 


2 

o>0 


3 

0 


/? 


(77) 


where  co0  is  the  orbital  rate,  and  making  the  small  angle  approximations,  the  gravity 
gradient  moment  equations  become 


Mg  =  3w0 


(/» 

-  y</> 

(/« 

-  wo 

0 

(78) 
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Equation(78)  clearly  shows  that  by  increasing  the  moments  of  inertia  for  /,  and  Iy 
with  respect  to  /2,  the  gravity  gradient  moments  for  roll  and  pitch  will  increase.  By  ex- 
tending a  point  mass  along  the  z  axis,  a  sizeable  increase  in  the  moments  of  inertia  along 
the  x  and  y  axes  is  realized. 

The  magnitudes  of  the  gravity  gradient  torques  generated  at  different  altitudes  is 
given  by  Figure  15. 
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Figure  15.      Gravity  Gradient  torque  Magnitudes  vs  Altitude  (Ref.  3,  p.  58] 

E.    DESIGN  FOR  STABILITY 

1.     Stability  and  Frequency  Analysis 

Combining  the  rigid  body  dynamic  equations  (61-63)  with  the  gravity  gradient 
equation  (78)  yields  the  undamped  linear  dynamic  equations  for  roll,  pitch,  and  yaw: 


Roll     Ix4>  +  4a>5(/  -  Q4>  -  oj0{Ix  +  I2-  IM  =  0 


(79) 


Pilch     LB  +  3u0{Ix  -  L)0  =  0 


(80) 
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Yaw     l2i>  +  a>5(/  -  Ix)\j/  +  to0(Ix  +  /,  -  I)4>  =  0 


(SI! 


Using  Reference  9  notations  for  simplicity,  the  following  variables  are  defined: 


A  = 


B  = 


I, 

— 

h 

h 

h 

— 

h 

(82) 


(83) 


C 


Jy-h 


(84) 


Substituting  eqs(82-S4)  into  eqs(79-81)  yields  the  following  second  order  differential 
equations 


<j>  +  4Aoj20<t)  -  (1  -  A)q)q\Ij  =  0 


(85) 


e  +  iBcoie  =  o 


(86) 


(A  +  CcoqiA  +  U  -C)4>  =  0 


(87) 


Taking  the  Laplace  transform  of  eqs( 85-87)  allows  for  frequency  analysis.  The  results 
are  put  into  matrix  form: 


(88) 


The  determinant  of  the  matrix  is  the  system's  characteristic  equation.  Using  the  middle 
row  to  simplify  calculations,  the  characteristic  equation  is  found  to  be 


S2  +  4Aa}20            0 

-(l-A)o)0S 

V 

"o 

0            S2  +  3Ba20 

0 

e 

= 

0 

{\-C)co20S           0 

S"  +  Cco0 

A 

_0 

(S4  +  (1  +  3  A  +  CA)(o20S2  +  4AC(Oq)(S2  +  3Bco20)  =  0 


(89) 


"Solution  of  the  first  term  on  the  left  yields  the  natural  frequencies  of  roll  and  yaw, 
repectively;  the  second  term  on  the  left  defines  the  natural  frequency  of  the  pitch  axis." 
[Ref.  9,  p.  9]  Solving  for  the  first  term  on  the  left  (roll/yaw  stability)  by  the  quadratic 
formula  and  dividing  through  by  col  yields 
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S?oii       Sl«w       ~(1  +  1A  +  AC)±  J(\  +  3A  +  AQ2-\6AC 


(90) 


For  roll 'yaw  stability  the  roots  of  eq(90)  must  be  real  and  negative;  therefore  the  fol- 
lowing stability  conditions  exist  [Ref.  4  p.  611] 


\  +  3A  +  AC>4jAC  (91) 

AC>0  (92) 

Solving  for  the  second  term  of  Eq(89)  and  dividing  through  by  cog  yields 

co0 

The  condition  for  stability  for  the  pitch  axis  is 

A>  C  (94) 

Converting  A.  B,  and  C  back  into  their  respective  moment  of  inertia  relationships  de- 
fined by  eqs(82-84).  combined  with  the  conditions  for  stability,  yields  two  possible  ori- 
entations: Iy>  Ix>  I2  and  I„>  I2>  Iy.  Figure  16  on  page  26  shows  a  summary  of  the 
stability  regions  based  on  moment  of  inertia  ratios.  [Ref.  10  p.  204] 
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Figure  16.      Gravity  Gradient  Stability  Regions 


26 


III.     ENVIRONMENTAL  DISTURBANCES  AND  RESTORING  TORQUES 

A.     ENVIRONMENTAL  DISTURBANCES/TORQUES 

The  are  four  major  environmental  disturbances  that  affect  an  orbiting  spacecraft: 
aerodynamic,  magnetic,  gravity  gradient,  and  solar  radiation  torque.  Figure  17  lists 
these  disturbances  and  their  relative  regions  of  dominance. 


SOURCE 

DEPENDENCE  ON  DISTANCE 

REGION  Or  SPACE  WHERE 

FROM  EARTH 

DOMINANT* 

AEROOYNAMIC 

or 
1 
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f. 

Figure   17.      Environmental  Disturbance  Torques  [Ref.  4,  p.  17] 

1.     Aerodynamic  disturbance 
a.     Drag 

"For  satellites  near  the  earth  the  principle  non-gravitational  force  is  aero- 
dynamic drag.  Aerodynamic  drag  is  a  retarding  force  due  to  atmospheric  friction  and  is 
in  the  direction  opposite  the  space  vehicles  velocity  vector."  [Ref.  4,  p.  63] 

The  drag  force  is  a  function  of  vehicle  velocity,  air  density,  and  surface  area 
of  the  satellite.  The  drag  force  is  given  by  the  following  equation: 


d  = 


pv  cja 


(95) 


where  d  is  the  drag  force,  p  is  the  atmospheric  density,  v  is  the  velocity  of  the  satellite 
along  its  orbital  path,  cd  is  the  drag  coefficient,  and  a  is  the  area  of  the  satellite  over 
which  the  drag  acts. 

b.     Aerodynamic  torque 

The  collisions  of  air  molecules  of  the  upper  atmosphere  with  the  satellite 
surface  produces  a  torque  about  its  center  of  mass.    Aerodynamic  disturbances  are  a 
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function  of  the  satellite's  altitude,  velocity,  and  symmetry.  A  satellite  with  a  relatively 
low  earth  orbit  may  be  significantly  affected  by  aerodynamic  torques. 

Shown  in  Figure  18  is  a  cylindrical  satellite  with  its  center  of  mass  (cm)  and 
center  of  pressure  (cp).  The  offset  between  the  center  of  mass  (cm)  and  the  center  of 
pressure  (cp)  is  denoted  by  Lx,  Iv,  and  L,,  respectively.  The  simplified  aerodynamic  dis- 
turbance torques  from  Ref  8  may  be  computed  by  the  following  expression: 


Ta  -  PJaAa  sin  a 


(96) 


where  Ta  is  the  torque  due  to  aerodynamic  pressure  (ft-lbs),  Pa  is  the  aerodynamic  pres- 
sure (lbs///2),  /„  is  the  distance  between  the  center  of  mass  and  the  center  of  pressure  (ft), 
Aa  is  the  exposed  surface  area  {ft7)  and  a  is  the  angle  of  attack  (radians).  For  practical 
purposes,  the  aerodynamic  torque  in  the  x  direction  is  zero  (sin  a  =  0  when  a»0).  The 
conversion  factor  from  ft-lbs  to  N-m  is  approximately  1.356.  Figure  19  on  page  29 
shows  aerodynamic  pressure  as  a  function  of  altitude. 


Figure  18.      Aerodynamic  effects  for  cylindrical  satellite  (Ref.  3,  p.  58] 
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\                               NOTE:      FROM  ATMOSPHERIC  DENSITY 

\                                                BASED  UPON  EXTENDED 
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Figure  19.      Aerodynamic  pressure  as  a  function  of  altitude  [Ref.  8,  p.  456] 

The  magnitude  of  the  pitch  and  yaw  aerodynamic  disturbances  for  a  one 
inch  cm-cp  offset  are  shown  in  Figure  20  on  page  30  for  various  altitudes.  The  magni- 
tude of  these  torques  decrease  significantly  with  increases  in  altitude.  "Proper  control  of 
the  satellites  inertias  to  enhance  gravity  gradient  torques  can  make  aerodynamic  dis- 
turbances essentially  insignificant  for  orbital  altitudes  above  300  nautical  miles  (480km). 
On  the  other  hand,  in  the  100  to  200  nautical  mile  (160-320km)  altitudes  aerodynamic 
disturbances  are  3  to  4  orders  of  magnitude  greater  than  the  gravity  gradient  torque  from 
a  1  slug-fr2  inertia  difference  and  1  degree  attitude  error.  In  the  lower  regions,  very  large 
inertia  differences  are  required  for  accuracy.  [This  leads  to  design  considerations]  where 
proper  configuration  can  make  these  aerodynamic  torques  work  as  restoring  torques 
instead  of  disturbance  torques"  [Ref.  3,  p.  57]. 
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Figure  20.      Aerodynamic  torque  for  a  1  inch  offset  as  a  function  of  altitude. 

2.      Rigid  Body  Solar  Pressure  Torque 

"The  major  factors  determining  the  radiation  torque  on  a  spacecraft  are  the  in- 
tensity and  spectral  distribution  of  the  incident  radiation,  the  geometry  of  the  surface 
and  its  optical  properties,  and  the  orientation  of  the  sun  vector  relative  to  the  spacecraft" 
[Ref.  4,  p.  570]. 

Solar  pressure  torques  are  the  disturbances  produced  by  solar  radiation  pres- 
sure and  is  a  function  of  the  offset  of  the  center  of  pressure  (cp)  and  the  center  of  mass 
(cm).  Referring  to  Figure  21  on  page  31,  the  center  of  mass-center  of  pressure  offsets 
are  designated  as  Lx.  Ly,  and  Lz,  respectively.  The  torques  generated  due  to  solar  pres- 
sure are  given  by  [Ref.  8,  p.  454]  as: 


T5  =  PSL5AS  cos  c 


Ts  =  2PSLSAS  cos  c 


(97) 
(98) 


where  T,  =  torque  due  to  solar  radiation  (ft-lbs),  P,  =  radiation  pressure  (lb/fr2)  and  is 
«  constant  at  9.65  x  10"  \b,'ft2  for  an  earth  orbiting  vehicle,  L,  =  center  of  mass  center 
of  pressure  offset  (ft),  As  =  surface  area  of  satellite  normal  to  sun(//2),  and  c  =  angle 
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of  incidence  (degrees).     Note,  equation  (97)  is  for  an  absorbent  satellite  body  and 
equation  (98)  is  for  a  reflective  body. 
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Rigid  body  solar  pressure  torques  [Ref.  3,  p.  61] 


B.     RESTORING  TORQUES 

As  seen  from  chapter  2,  the  restoring  torques  of  the  undamped  rigid  body  combined 
with  the  gravity  gradient  effects  are  considerably  smaller  for  the  yaw  axis  than  for  either 
the  pitch  or  roll  axes.  This  is  due  to  the  necessity  of  large  moments  of  inertia  along  the 
pitch  and  roll  axes,  as  compared  to  the  yaw  axis,  for  gravity  gradient  stabilization  to 
occur.  The  lack  of  adequate  yaw  restoring  torque  requires  that  additional  yaw  restoring 
torque  be  provided  for  accurate  3-axis  stabilization.  A  discussion  of  yaw  restoring  torque 
methods  follows. 

1.     Thrusters 

Thrusters  are  an  effective  torque  generating  device;  however,  for  the  purpose 
of  this  thesis,  thrusters  will  not  be  considered  for  attitude  control.  Thrusters  are  neces- 
sary for  orbit  maintenance  and  orbital  transfer.  Due  to  the  physical  size  constraints  im- 
posed by  the  GAS  cannister,  ORION  is  severely  limited  in  its  fuel  storage  capacity  to 
approximately  71  lbs  of  hydrazine  fuel  [Ref  1].  An  assumption  that  the  launch  platform 
will  place  ORION  in  its  final  orbit  is  optimistic.  Figure  22  on  page  32  shows  the  fuel 
required  to  go  from  one  orbital  altitude  to  another.  Saving  the  onboard  fuel  for  orbit 
maintenance  and  transfer  and  not  for  attitude  control  allows  more  fuel  to  be  devoted 
towards  supporting  the  mission  either  by  allowing  greater  orbital  flexibility  and/or  a 
longer  lifetime. 
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Figure  22.      Fuel  required  for  orbit  transfer  based  on  a  2501b  satellite.  [Ref  2,  p.  23] 

2.     Momentum  Exchange  Devices 

The  principle  behind  momentum  exchange  devices  is  that  by  changing  the  an- 
gular momentum  of  the  control  device,  the  angular  momentum  of  the  vehicle  will  be 
changed  an  equal  and  opposite  amount.  Recall  that  the  rate  of  change  of  angular  mo- 
mentum is  equal  to  the  torque  that  is  generated.  The  two  basic  types  of  momentum  ex- 
change devices  are  reaction  wheels  and  control  moment  gyros  (CMG).  The  basic 
difference  between  the  two  is  the  method  by  which  the  change  in  angular  momentum  is 
accomplished.  The  reaction  wheel  axis  is  fixed  (generally  along  one  of  the  axes  of  the 
vehicle)  so  that  the  change  in  angular  momentum  is  accomplished  by  varying  the  speed 
of  rotation  of  the  flywheel.  The  CMG  creates  the  change  in  angular  momentum  by  tilt- 
ing the  contant  speed  flywheel  with  respect  to  the  vehicle. 
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a.     Reaction  wheels 

"An  axis  may  be  controlled  by  varying  the  speed  of  the  reaction  wheel  in 
response  to  an  attitude  error  [Ref.  6,  p.  149].  The  momentum  of  a  system  with  a  reaction 
wheel  is 


77  =  Hw  +  Hc 


(99) 


where  77c  is  the  angular  momentum  of  the  system  and 


Hw  = 


77, 


wy 


Hu 


is  the  angluar  momentum  of  the  reaction  wheel.  Aligning  the  reaction  wheels  along  each 
axis,  and  substituting  eq(99)  into  eq(48)  and  adding  the  gravity  gradient  moments  of 
eq(7S)  yields 

ljj>  +  Aco20(Iy  -  Iz)4>  -  co0{Ix  +  I2  -  Iy)>P  +  (6-  co0)HW2  -  ty  +  aj0<t))Hwy  +  Hwx  =  JY, 

Iy6  +  30)5(7,  -  I2)6  +  (  -  4>  +  coQ^)HW2  +  (<A  +  coQ4>)Hwx  +  Hwy  =  Yjy  (101) 

I2'\jt  +  coq(/v  -  7x)i//  +  Q)0(IX  +  I2  -  Iy)4>  +  (</>-  co0\J/)  +  (-</>  +  o)0\j/)Hwx  +  Hw.  =  l_Jz 

"Because  Hwx,Hwy,HW2,  and  co0  are  small,  the  coupling  terms  are  small.  If  the  coupling 
terms  are  neglected,  the  equations  of  motion  about  the  roll,  pitch,  and  yaw  axes  become 
independent,  and  hence,  they  can  be  controlled  independently"  [Ref.  6,  pp.  149-150]. 
Dropping  the  coupled  terms,  the  motion  equations  become 


Ix4>  +  A(ol{Iy  -  I2)4>  -  co0(7x  +  I2-  ly)4>  +  Hwx  =  Y,TX 


(102) 


Iyd  +  30)2(7,  -  I2)d  +  Hwy  -  Yjy 


(103) 


72«A  +  u)0(Iy  -Ix)ii>  +  co0(Ix  +  IZ-  L)4>  +  HW2  =  >  T2 


(104) 
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Allowing  the  momentum  of  the  reaction  wheel  to  be  proportional  to  the  attitude  error, 
then  by  differentiation,  the  reaction  wheel  control  torques  are  equal  to 


Txm  =  -  Kx(rxj>  +  (f>) 
Tyrw  =  -  ky(ry6  +  6) 

T2rw  =  ~  ^(vA  +  <A) 


(105) 


(106) 


(107) 


and  the  control  scheme  would  be  as  shown  in  Figure  23. 
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Figure  23.      3  axis  reaction  wheel  control  scheme  [Ref.  6,  p.  149] 

The  torque  genterated  by  the  reaction  wheel  is  a  function  of  the  wheel's  angular  velocity 
and  its  moment  of  inertia.  This  means  that  a  small  lightweight  wheel  at  a  high  rate  of 
speed  generates  the  same  amount  of  torque  that  a  large  wheel  at  a  low  rate  of  speed. 
Since  ORION  is  designed  for  a  relatively  short  lifespan,  wear  and  tear  on  spinning  parts 
is  not  an  overriding  concern.  It  should  be  noted  that  continuous  disturbance  torques 
that  are  not  cyclical  in  nature  will  eventually  saturate  the  reaction  wheel,  whereby,  the 
reaction  wheel  rotational  speed  will  reach  a  maximum.  Thrusters  or  some  other  method 
of  momentum  removal  must  be  used  to  reduce  wheel  speed.  Disturbance  torques  that 
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are  cyclical  in  nature  arc  easily  offset  by  the  reaction  wheel  as  long  as  the  disturbance 

torque  is  such  that  it  does  not  saturate  the  reaction  wheel  speed. 

3.     CMG(s) 

"A  control  moment  gyro's  (either  single  or  double  gimbaled)  angular  mo- 
mentum is  due  to  the  rotor  which  is  spinning  about  the  spin  axis  with  a  constant 
angular  rate.  Because  the  spin  axis  is  gimbaled,  a  commanded  gimbal  rotation  causes 
the  direction  of  the  angular  momentum  vector  to  change,  thus  creating  a  control 
torque  parallel  to  the  output  axis.  I  he  magnitude  of  this  torque  depends  on  the 
speed  of  the  rotor  and  the  gimbal  rotation  rate."   [Ref.  4,  p.  200] 

A  control  moment  gyro  scheme  is  shown  in  Figure  24. 


Figure  24.      T>\in  gyro  controller  with  one  drive  motor  [Ref.  8,  p.  415] 

Because  of  their  expense  and  weight,  CMG's  are  only  used  on  large  spacecraft 
[Ref.  4].    For  this  reason,  control  moment  gyros  will  not  be  addressed. 

4.     Magnetic  torquers 

Magnetic  torquers  may  be  used  to  dampen  the  kinetic  energy  of  the  libration 
motion  as  well  as  correcting  attitude  errors  in  any  direction.  The  magnetic  torquers 
consist  of  three  wire  coil  sets  arranged  perpendicularly  to  each  other.  This  configuration 
allows  magnetic  torque  to  be  generated  along  any  axis  or  in  any  direction  as  a  result  of 
adding  two  or  three  different  orthogonal  torques.  As  shown  in  Figure  25  on  page  36, 
when  a  current  is  applied  around  the  loop  of  a  coil,  it  produces  a  magnetic  dipole  that 
is  normal  to  the  plane  of  the  coil  and  a  magnitude  which  is  proportional  to  the  coil's 
enclosed  area  and  ampere-turns. 
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Figure  25.      Magnetic  Moment  Due  to  a  Current  Loop  [Ref.  4,  p.  204] 


The  magnetic  moment  M  is  defined  as 


M  =  NIAn 


(108) 


where  N  is  the  the  number  of  turns  around  the  coil,  I  is  the  current  in  amperes,  A  is  the 
area  enclosed  by  the  loop,  and  h  is  the  unit  vector  normal  to  the  plane  of  the  loop.  Be- 
cause of  the  relationship  between  N,  I,  and  A,  many  different  combinations  are  available 
to  produce  the  desired  magnitude  of  the  magnetic  moment.  Figure  26  on  page  37  shows 
three  existing  spacecraft  with  various  combinations  of  number  of  turns  and  current 
drawn.    The  torque  generated  by  the  magnetic  torquer  is 


Tmas  =  MxB 


(109) 


where  M  is  the  magnetic  dipole  and  B  is  the  earth's  magnetic  field.  The  earth's  magnetic 
field  for  an  equatorial  orbit  can  be  approximately  modeled  as  a  tilted  dipole  located  at 
78.3°  north  latitude  and  69°  west  longitude  as  shown  in  Figure  27  on  page  37. 
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Figure  26.      Some  existing  magnetic  moment  schemes  [Ref.  3,  p.  205] 
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Figure  27.      Earth's  magnetic  field  modeled  as  a  tilted  dipole  [Ref.  6,  p.  148] 

The  strength  of  the  earth's  magnetic  field  for  a  circular  equatorial  orbit  in  component 
form  can  be  approximated  by 
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Bx  = f-sin(11.7)sin(281  -  lams)  (110) 

R 

Mg 

By  = f-cos(11.7)  (111) 

R3 

M. 
B2  =  — f-sin(11.7)cos(281  -  lams)  (112) 

R 

where  lami  is  measured  in  degrees  of  right  ascension  from  Grenich,  Mt  is  the  strength 
at  the  earth's  center  which  is  approximately  8  x  101S  Wb-m,  and  R  is  the  distance  in 
meters  from  the  earth's  center  to  the  orbital  altitude.  With  respect  to  the  satellite,  the 
earth's  magnetic  field  is  changing  continuously  as  the  vehicle  moves  along  its  orbital 
path.   Putting  into  component  form  both  M  and  B,  eq(109)  becomes 

Tmag_x  =  M2By-MyB2  (113) 

Tmag_y  =  MxB2-M2Bx  (114) 

Tmag_2  =  MyBx-MxBy  (115) 

It  should  be  noted  that  the  earth's  magnetic  field  is  not  uniform  and  varies  in 
direction  based  on  the  satellite's  position.  For  accurate  attitude  control,  the  vehicle  must 
be  able  to  exactly  know  both  its  position  and  the  earth's  magnetic  field  at  that  position. 
For  low  earth  orbits,  the  magnetic  field  is  well  mapped  and  may  be  easily  programmed 
into  the  attitude  control  computer.  For  the  purpose  of  this  thesis,  magnetic  field  ap- 
proximations will  suffice  to  demonstrate  magnetic  torquing  as  a  viable  attitude  control 
capability. 
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IV.      CHAPTER  4 

A.      GRAVITY  GRADIENT  STABILIZATION  RESULTS 

1.  Attitude  Control 

There  are  three  general  attitude  control  phases  with  respect  to  gravity  gradient 
systems:  gravity  gradient  capture,  transition  period,  and  steady  state.  Gravity  gradient 
from  a  tumbling  mode  occurs  when  angular  velocity  slows  to  one  rotation  per  orbit  or 
less,  allowing  the  minimum  moment  of  inertia  axis  to  align  with  the  local  vertical.  The 
transition  period  is  the  time  from  gravity  gradient  capture  until  a  steady  state  has  been 
achieved.  This  period  of  oscillations  will  vary  depending  on  the  dampening  scheme  used 
to  reduce  the  amplitude  of  oscillations.  Predominant  throughout  the  satellite's  lifetime, 
the  steady  state  phase  is  the  period  in  which  the  satellite's  errors  are  confined  within 
operationally  acceptable  amplitudes.  Maintenance  of  the  desired  steady  state  requires 
minimization  of  external  disturbances  and  control  of  those  which  cannot  be  eliminated. 

The  following  group  of  simulations  assumes  that  a  steady  state  phase  has  been 
achieved  and  that  the  small  angle  approximations  are  in  effect.  For  each  of  the  simu- 
lations, the  following  initial  conditions  are  in  effect:  <b,  8,  and  \\i  position  error  =  5°,  and 
<b,  6  ,  and  4/  =  0. 

2.  Various  results  based  on  different  configurations 

a.  Effects  of  extending  booms 

For  Figure  28  on  page  40,  a  10  kg  point  mass  is  used  and  the  mass  of  the 
boom  is  assumed  to  be  negligible.  Extension  of  the  boom  increases  the  moment  of  in- 
ertia in  the  x  and  y  directions.  This  increases  the  restoring  torque  along  the  roll  and  pitch 
axes  although  it  does  little  to  dampen  the  libration  motion. 

b.  Gravity  gradient  effect  alone 

The  following  figures  show  the  gravity  gradient  effect  on  a  satellite  at  an 
altitude  of  1000  km  having  moments  of  inertias  as  shown.  As  shown  by  Figure  29  on 
page  41  and  Figure  30  on  page  42,  stabilization  along  the  roll  and  pitch  axes  would  take 
a  prohibitively  long  time  while  the  yaw  direction  would  not  stabilize.  It  is  clear  that 
gravity  gradient  stabilization  alone  is  incapable  of  acceptable  three  axis  stabilization. 
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MOMENT  OF  INERTIA  VS  DOOM  LENGTH 
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Figure  28. 


Effect  of  extending  boom  on  moments  of  inertia 
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Figure  29.      Gravity  gradient  effect  only  Ix=  15,  Iy=  18,  Iz=  3 
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Figure  30.      Gravity  gradient  effect  only  Ix  =  95,  Iy  =  98,  Iz=  3 


c.      Effects  of  Reaction  Wheel(s) 


One  method  of  augmenting  gravity  gradient  stabilization  is  through  the  use 
of  reaction  wheels.  Teldrix  reaction  wheel  classes  and  corresponding  technical  data  are 
presented  in  Figure  31  on  page  43. 
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Wheel  diameter 

r.m 

Nins 

20 
1.8. .6. 5 

26 
5.0. .20 

35 

50 
50.. .300 

60') 
200...  1000 

Angular  momentum  tango 

14. .80 

Max.  reaction  torque 

Nm 

0.2 

0.2 

0.2 

0.3 

0.3-0.6 

Speed'") 

min  ' 

6000 

6000 

6000 

6000 

6000 

Loss  torque  al  max.  speed") 

Nm 

<  0.012 

<  0.013 

<  0.015 

<  0.022 

<0.07 

Power  consumption: 

-  steady  stale  (depending 
on  speed) 

W 

2.7 

2. .8 

2... 10 

3... 15 

10. .50 

—  max.  power  rating 

W 

<60 

<80 

<  too 

<  150 

<500 

Dimensions: 

—  diameter  A 

mm 

203 

260 

350 

500 

600 

-  height  B 

mm 

75 

85 

120 

150 

180 

Weight 

^ 

2. 7.  .3.-1 

3.5.6.0 

5.0...B.0 

7.5...12 

20. .37 

Environmental  conditions: 

-  operating  temperature 

-  vibration  (sinusoidal) 

-  vibration  (random) 

-  linear  acceleration 

suitable  lor  satellites  compatible  with  launchers 
such  as  ARIANE  or  Space  Shuttle 

Figure  31.      Teldiix  Reaction  Wheel  Technical  Data 

Of  the  five  diameter  classes,  the  20  cm  diameter  wheel  is  the  most  appropriate  of  the 
classes  listed  for  ORION.  Reaction  wheels  also  come  in  diameters  smaller  than  20  cm 
diameter.  Bendix  Aerospace  SP  reaction  wheel  of  Figure  32  on  page  44  has  a  diameter 
of  6  inches. 

As  illustrated  in  Figure  31  and  Figure  32  on  page  44,  power  consumption 
values  differ  with  dilfercnt  reaction  wheels.  For  the  20  cm  Teldrix  reaction  wheel,  steady 
state  requires  2-7  watts  with  a  maximum  less  than  60  watts.  The  Bendix  SP  reaction 
wheel  requires  2.5  watts  at  steady  state  with  a  maximum  power  consumption  of  6.5 
watts.  Due  to  weight  and  power  consumption  values,  the  Bendix  reaction  wheel  appears 
to  be  the  superior  choice  of  the  two. 

The  following  figures  were  generated  by  setting  the  reaction  wheel  torques 
equal  to  eqs(107-109).  The  initial  conditions  are  0  ,  6,  and  \j/  position  errors  =  5°,  and 
<f>  =  6  =  \\j  =  0.  For  simulation  purposes,  the  maximum  reaction  torque  value  of  .2  N-m 
from  Figure  31  will  be  in  effect. 

Figures  33-35  have  moment  of  inertias  /,  =  15,  /,=  18,  and  /,  =  3,  with 
K,  =  Ky  =  2.0,  K2  =  1.0,  and  t,  =  xy  =  x2  =  1.0.  Figure  33  on  page  45  shows  <fi  (position 
error),  <f>  (velocity),  and  the  restoring  torque  (Trw6)  generated  by  the  reaction  wheel  with 
respect  to  time.  Figures  34  and  35  show  6,  6  ,  7^,  and  \{/,  ty  ,  Trwt,  respectively.  With 
the  relatively  small  moments  of  inertias,  stabilization  of  all  axes  within  the  ±  1°  criteria 
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occurs  within  40  seconds.  Increasing  the  moments  of  inertias  to  Ix  =  95,  /,,  =  98  ,  and 
I2  =  3  while  holding  A'„  Ky  at  2.0,  A',  at  1.0,  and  t„  t,  ,  and  t2  constant  at  1.0,  yields 
Figures  36-38.  As  expected,  the  </>  and  0  axes  take  longer  to  stabilized  (approximately 
160  seconds)  than  the  above  case  with  the  smaller  moments  of  inertia.  The  \p  axis  of 
Figure  38  stabilizes  much  faster  than  either  the  </>  or  0  axes.  This  is  due  to  the  relatively 
small  moment  of  inertia  (/z)  compared  to  the  restoring  torque  generated  by  the  reaction 
wheel. 
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Figure  32.      BendLx  Aerospace  Reaction  Wheel  Technical  Data. 
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Figure  33.       Gravity  gradient  effect  with  3  RVV  (X  axis)  L\=15,  Iy=18,  Iz=3 

K,  =  2.0,t,=  1.0 
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Figure  34.       Gravity  gradient  effect  with  3  RVV  (Y  axis)  Ix=15,  Iy=  18,  Iz=3 

Ky=2.0  ,ry=  1.0 
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Figure  35.       Gravity  gradient  effect  with  3  RVV  (Z  axis)  Ix=  15,  ly  =  IS,  Iz=3 

a;  =  2.o,t2=  i.o 
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Figure  36.       Gravity  gradient  effect  with  3  RVV  (X  axis)  Ix=95,  Iy  =  98,  lz=3 

A',  =  2.0  ,  t,  =  1.0 
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Figure  37.       Gravity  gradient  effect  with  3  RW  (Y  axis)  Ix  =  95,  Iy  =  98,  Iz=3 
^  =  2.0,t,=  1.0 
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Figure  38.       Gravity  gradient  effect  with  3  RVV  (Z  axis)  Ix=95,  Iy  =  98,  Iz=3 

K,=  1.0  ,  t2=  1.0 
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Figures  39-41  have  moments  of  inertia  of  Ix  =  15,  Iy  =  18,  and  I,  =  3  with 
Kx  =  Ky  =  .2,  K2  =  .1,  and  tx  —  ry  =  r2  =  1.0.  The  reduction  of  Kx,  Ky  and  A.',  reduces  the 
reaction  wheel  torque  generated. 

Figure  39  shows  ^(position),  0( velocity),  and  reaction  wheel  torque  (7^) 
with  respect  to  time.  Figures  40  and  41  show  6,  6,  TrwS  ,  and  \}/,  <//,  Tw  ,  respectively. 
Stabilization,  within  the  ±  1°  criteria,  occurs  for  all  axes  in  under  200  seconds.  Com- 
parison of  the  magnitudes  of  the  torque  values  of  Figures  33-35  and  Figures  39-41  for 
the  same  moment  of  inertia  values  shows  a  decrease  of  roughly  a  factor  of  five. 

Increasing  the  moments  of  inertia  to  Ix  =  95,  Iy  =  98  and  I2  =  3,  while  hold- 
ing Kx  =  Ky  =  .2,  K2  =  .1,  and  t„  =  t,  =  t.  =  1.0,  produces  Figures  42-44.  This  last  group 
of  reaction  wheel  figures  shows  that  s:650  seconds  are  required  to  dampen  the  system 
to  the  ±1°  stabilization  criteria.  The  tradeoffs  between  the  different  reaction  wheel 
combinations  are  essentially  time  response  against  torque/power.  The  faster  the  desired 
response,  the  greater  the  magnitude  of  the  restoring  torque,  and  ultimately,  the  more 
power  that  will  be  consumed  in  achieving  this  response. 
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Figure  39.       Gravity  gradient  effect  with  3  RW  (X  axis)  Ix=  15,  Iy=  18,  Iz=3 


K,  =  .2  ,  t,  =  1.0 
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Figure  40.       Gravity  gradient  effect  with  3  RVV  (Y  axis)  Ix=  15,  iy-  18,  Iz=3 

Ky=  .2  ,t,=  1.0 


53 


T 


RW^ 


Figure  41.       Gravity  gradient  effect  with  3  RVV  (Z  axis)  Ix=15,  ly  =  18,  Iz=3 
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Figure  42.       Gravity  gradient  effect  with  3  RW  (X  axis)  Ix  =  95,  Iy  =  98,  Iz=3 
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Figure  43.       Gravity  gradient  effect  with  3  (Y  axis)  RW  Ix=95,  Iy  =  98,  Iz=3 

a;  =  .2 ,  t,=  l.o 


56 


s 


r- 


O    n 


•f  IM  t4t  HI  IH  «H«  5M  «4t  T7t  Ml 

rtrt  imci 


r. 


R\V> 


IM  m  4ii  •»  sw  m  m 

Tire  iseci 


Figure  44.       Gravity  gradient  effect  with  3  RVV  (Z  axis)  Ix  =  95,  Iy  =  98,  Iz=3 
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d.      Effects  of  Magnetic  torquers 

For  the  magnetic  torquing  simulations,  the  initial  conditions  are  <f>,  6,  and 
\jj  position  errors  =  5°  ,  and  <fi.  6,  and  \J/  =  0.  The  initial  restoring  torque  is  provided 
by  the  gravity  gradient  effect  combined  with  the  rigid  body  dynamics.  For  the  following 
simulations,  one  of  the  two  conditions  must  exist  prior  to  energizing  the  magnetic  coils 
to  generate  magnetic  torque: 

1)  position  greater  than  .5°  and  velocity  greater  than  0,  or 

2)  position  less  than  —.5°  and  velocity  less  than  0.  For  the  following  figures, 
the  magnetic  moment  varies  from  9.0427  to  .5  Wb-m.  At  the  maximum  value  of 
9.0427W-m,  N=400  turns.  I  =  .5  amp,  and  A  =  nr2  where  r  is  the  radius  at  .12m.  Subse- 
quent values  are  obtained  by  reducing  the  current  drawn. 

The  magnetic  Field  in  component  form  for  an  equatorial  orbit  can  be  ap- 
proximated by  eqs(  1 10-1 12).  Since  the  earth's  magnetic  field  is  highly  dependent  on  the 
satellite's  position,  ORION  must  have  the  capability  to  determine  its  position  in  order 
to  effectively  utilize  magnetic  torquing.  For  simulation  purposes,  the  satellite  is  assumed 
to  be  directly  over  Greenwich  at  time  t  =  0.  The  magnetic  coils  are  turned  on  and  off  at 
the  specified  strength  as  needed.  Varying  the  strengths  of  the  magnetic  moments  gener- 
ated would  enhance  performance,  but  is  beyond  the  scope  of  this  thesis.  The  assumption 
that  magnetic  torque  may  be  generated  in  any  direction  due  to  vector  addition  is  not 
always  valid:  magnetic  torque  may  not  be  generated  in  the  direction  of  the  magnetic 
field. 

Fieures  45-47  utilized  magnetic  moments  of  m_  =  mv  =  m  =  9.0427  with 
moments  of  inertias  Ix=  15,  Iy=  18,  and  l7  =  3.  For  Figure  45  on  page  60,  the  gravity 
gradient  restoring  torque  takes  approximately  800  seconds  before  meeting  a  condition 
necessary  for  energizing  the  magnetic  coils.  Due  to  the  small  moments  of  inertia,  the  4> 
restoring  torque  is  needed  for  only  a  short  time.  Figures  46  and  47  show  the  6  and  \p 
axes  essentially  stabilized  within  the  ±  1°  criteria  after  the  the  initial  magnetic  torque 
pulse.  For  the  6  axis,  this  pulse  occurs  at  the  1200  seconds  and  lasts  for  approximately 
300  seconds.  The  d  axis  has  a  lower  frequency  response  and,  therefore,  takes  longer  to 
to  meet  the  condition  necessary  to  energize  the  magnetic  coils.  Due  to  its  relatively  small 
moment  of  inertia,  the  \j/  axis  chatters  slightly  before  stabilizing.  A  smaller  magnitude 
of  restoring  torque  eliminates  this  problem. 

Increasing  the  moments  of  inertia  to  Ix  =  95,  Iy  =  98,  and  I,  =  3  while  leav- 
ing the  magnetic  moments  mx,  my,  and  mz  constant  at  9.0427wb-m,  produces  Figures 
48-50.  As  expected,  more  restoring  torque  is  needed  to  stabilized  the  <j>  and  6  axes  forc- 
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ing  the  magnetic  coils  to  remain  energized  for  a  longer  period  of  time.  The  time  needed 
for  all  axes  to  be  within  the  ±  1°  criteria  is  limited  by  the  6  axis.  Stabilization  occurs 
within  2800  seconds  while  for  the  case  with  the  smaller  moments  of  inertias,  stabilization 
occurs  at  approximately  2100  seconds. 
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Figure  45.        Gravity    gradient    with    magnetic    torquing    (X    axis):     Magnetic 
moment  =9.0427W-m  Ix=  15,  Iy=  18,  Iz=  3. 
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Figure  46.        Gravity    gradient    with    magnetic    torquing    (Y    axis):     iMagnetic 
moment  =  9.0427W-m  Ix=  15,  Iy=  18,  Iz=3. 
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Figure  47.        Gravity    gradient     with    magnetic    torquing    (Z    axis):     Magnetic 
moment  =  9.0427W-m  Ix=  15,  Iy=  18,  Iz=3. 
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Figure  48.        Gravity    gradient    with    magnetic    torquing    (X    axis):     iMagnetic 
moment  =  9.0427W-m  Ix  =  95,  Iy  =  98,  Iz=3. 
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Figure  49.        Gravity    gradient     with    magnetic    torquing    (Y    axis):     Magnetic 
moment  =9.0427W-mIx=  95,  Iy=98,  Iz=3. 
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Figure  50.        Gravity    gradient    with    magnetic    torquing    (Z    axis):     Magnetic 
moment  =9.0427W-m  Ix=95,  Iy  =  98,  Iz=3. 
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Reducing  the  strength  of  the  magnetic  moments  reduces  the  torque  that  is 
generated  and.  thereby,  slows  the  response  time  to  stabilization.  For  Figures  51-56,  the 
magnetic  moments  mx  =  my  =  mz  =  1.0.  Figures  51-53  have  moment  of  inertias  /,  =  15. 
ly  =  IS  ,  and  /,  =  3.  Stabilization  for  the  6  axis  occurs  within  approximately  3600  sec- 
onds. Increasing  the  moments  of  inertias  to  Ix  =  95,  Iy  =  98,  and  Iz  =  3  produces  Figures 
54-56.  Stabilization  within  the  +  1°  criteria  for  the  0  axis  occurs  at  approximately  11000 
seconds.  The  length  of  time  that  the  magnetic  torqers  are  energized  increases  as  the 
strength  of  the  moment  decreases.  Figures  57-62  show  the  effects  of  magnetic  torquing 
using  a  maximum  magnetic  moment  of  .5Wb-m  for  each  coil.  Figures  57-59  show  that 
for  low  moments  of  inertia  (/,  =  15,  [y  =  18,  and  72  =  3),  stabilization  occurs  within  the 
±  1°  within  4800  seconds.  Increasing  the  moments  of  inertia  to  Ix  =  95,  Iy  =  9S,  and 
I2  —  3  holding  the  maximum  magnetic  moment  =  .5\Vb-m,  produces  Figures  60-62.  In 
Figure  61,  the  B  axis  may  take  a  prohibitively  long  time. 


66 


* 


nillttCI   IS  1l«    IGUIIBtKI    IX-1S   It-tB   U-3.89 


4> 


i  «  Ti  M  I.I  i.l  i.i  I.I  M  II  >••  ••• 

lire  in'  SCO 

niiOIWBMCI   »5  UK    IKI1DMKI    IX-IS   IT-H   U-J.i 


T. 


mag—x 


1.1        I.I         I.I         I.I         I.I         I.I         I.I         I.I         1.1         1.1         I.I 

IIH   •!»•  SCI 

i«  v<;  i  nr  irr-]i*nMi 


Figure  51.        Gravity    gradient    with    magnetic    torquing    (X    axis):     iMagnetic 
moment=  l.OW-m  Ix=  15,  Iy=  18,  Iz=3. 
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Figure  52.         Gravity     gradient     with    magnetic    torquing    (Y    axis):     Magnetic 
moment  =  l.OW-m  Ix=  15,  Iy=  18,  Iz=3. 


68 


« » 

« < 

it 

J 

f.« 

V 

'     ,. 

t.t 

i 

\ 

I  • 

■1.1. 

•. 

•            •  '            ••            1.?            l'.«            J.»            I.I            t.l            I.I 

iit  tic  so 

J.I        I.I 

rsnotci  »s  HHt  icwuoom  ix-ts  n-»e  il-3.li 

!• 

• 

\ 

'/' 

■•I 
i« 

\_ 

I 

■HI 

I 

■HI 

c 

i 

»•«. 

i 

•        •••         t.i         i.i         i'.i         i.i         i.i         t.i         it 

lltt   lit'  SCI 

I.I        t.l 

rsioiwxretci  vs  tin  iccaitoKi  i«-is  n-u  ii-m 

»• 

it 

T 

mag. 

-2              i 

i 

T 

• 

> 

•l 

n                                                                                     , 

l.t          I.i          t.l          i.i      .    i.i          i.i          i.i          t.l          I.I 
III*   lit*  Ktl 

I.I 

I.I 

U  »S  lirt    IGU10IOKI 

Figure  53.        Gravity    gradient    with    magnetic    torquing    (Z    axis):     Magnetic 
moment  =  l.OW-m  Ix=  15,  Iy=  18,  Iz=  3. 
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Figure  54.         Gravity     gradient     with    magnetic    torquing    (X    axis):     Magnetic 
moment  =  l.OW-m  Ix=95,  Iy=98,  Iz=3. 
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Figure  57.        Gravity    gradient     with    magnetic    torquing    (X    axis):     Magnetic 
moment  =.50W-m  Ix=  15,  Iy=  18,  lz=3. 
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Figure  58.         Gravity     gradient     with    magnetic    torquing    (Y    axis):     Magnetic 
moment  =  .50\V-m  Ix=  15,  Iy=  18,  lz=3. 
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Figure  59.        Gravity    gradient    with    magnetic    torquing    (Z    axis):     Magnetic 
moment  =  .50\V-m  Ix=  15,  Iy=  18,  Iz=  3. 
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Figure  60.        Gravity    gradient    with    magnetic    torquing    (X    AXIS):     Magnetic 
moment  =. 50W-m  lx=  95,  Iy=98,  Iz=3. 
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Figure  61.        Gravity    gradient    with    magnetic    torquing    (Y    AXIS):     Magnetic 
moment  =.50W-mIx=  95,  Iy  =  98,  lz=3. 
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Figure  62.        Gravity    gradient    with    magnetic    torquing    (Z    AXIS):     Magnetic 
moment  =.50W-m  Ix=95,  Iy  =  98,  Iz=  3. 
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3.      Gravity  Gradient  capture  (from  tumbling  mode) 

Increasing  the  inertia  will  reduce  the  tumbling  rate  by  conforming  to  the  con- 
servation of  angular  momentum,  70o>0  =  /,&.),.  As  the  inertia  increases,  a  corresponding 
decrease  in  angular  rate  must  be  realized.  This  is  useful  in  obtaining  the  initial  gravity 
gradient  capture  from  a  tumbling  mode.  Figure  63  on  page  80  shows  that  for  a  slow 
tumbling  satellite,  the  change  in  the  moment  of  inertia  by  extending  the  boom  may  be 
sufficient  to  allow  gravity  gradient  capture.  The  satellite  is  initially  tumbling  at  a  rate 
of  5  rotations  per  orbit  about  the  pitch  axis.  The  extension  of  the  boom  increases  the 
moment  of  inertia  about  the  y  axis  by  a  factor  of  6;  therefore,  by  the  conservation  of 
angular  momentum,  the  angular  rate  must  decrease  by  a  factor  of  6.  This  reduces  the 
rotations  per  orbit  to  approximately  5/6  per  orbit  well  below  the  one  rotation  per  orbit 
rate  necessary  to  acheive  gravity  gradient  capture.  Figure  63a  shows  the  initial  tumbling 
rate  of  5  rotations  per  orbit.  Figure  63b  shows  the  corresponding  decrease  in  the  tum- 
bling rate  as  the  moment  of  inertia  is  increased.  Figure  63c  shows  the  increase  in  mo- 
ment of  inertia  as  the  boom  is  extended  with  respect  to  time. 

When  an  initial  tumbling  rate  is  higher  than  the  rate  that  can  be  captured  solely 
through  the  extension  of  the  boom,  magnetic  torquing  may  aid  in  gravity  gradient  cap- 
ture. By  creating  torque  in  the  direction  opposite  to  that  of  the  angular  velocity,  mag- 
netic torquing  can  be  used  to  slow  the  rate  of  the  tumbling  vehicle. 
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Figure  63.      Gravity  gradient  capture  by  extension  of  boom 
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CHAPTER  5 


A.     CONCLUSIONS 

1.  Gravity  gradient  stabilization  (only) 

Gravity  gradient  stabilization  alone  is  not  a  viable  attitude  control  method  for 
ORION*.  This  is  due  to  gravity  gradient  effect  providing  little  yaw  restoring  torque.  In 
addition,  gravity  gradient  stabilization  provides  little  dampening  efiect  on  the  roll  and 
pitch  axes.  When  augmented  by  other  control  schemes,  gravity  gradient  stabilization 
provides  sufficient  passive  restoring  torque  to  warrant  its  use. 

2.  Gravity  gradient  with  3  reaction  wheels 

a.  Advantages 

Gravity  gradient  stabilization  augmented  by  three  reaction  wheels,  each 
aligned  along  the  body's  axes,  provides  excellent  attitude  control  results  both  in  time 
response  and  accuracy. 

b.  Disadvantages 

Power  consumption  may  be  prohibitively  high  for  this  method.  At  steady 
state,  each  reaction  may  use  up  to  6  watts  while  drawing  considerably  more  during  the 
acceleration  or  deceleration  phase.  If  a  steady  state  drain  of  18  watts  (6  for  each  re- 
action wheel)  is  acceptable,  then  this  method  should  be  considered.  The  projected  end- 
of-life  power  budget  for  ORION  is  60  watts  [Ref.  1,  p.  132J.  Using  this  figure  as  an 
estimate,  a  steady  state  draw  of  18  watts  leaves  42  watts  available  for  mission  support. 
From  Figure  3  on  page  3,  an  average  of  34  watts  is  required  for  mission  support.  The 
size  and  weight  of  three  reaction  wheels  may  be  prohibitively  high  as  well.  Another 
disadvantage  is  that  thrusters  must  be  used  to  remove  momentum  when  the  reaction 
wheel(s)  become  saturated. 

3.  Control  moment  gyro  (CMG) 

a.  Advantages 

Although  not  discussed  in  detail,  control  moment  gyros  are  capable  of 
generating  large  restoring  torques. 

b.  Disadvantages 

The  weight  requirements  for  a  control  moment  gyro  are  prohibitively  high 
when  compared  to  reaction  wheels  when  relatively  small  torques  are  invoved.  As  a  small 


multi-purpose  satellite,  ORION  does  not  need  the  magnitude  of  restoring  torque  neces- 
sary to  warrant  the  use  of  control  moment  gyros. 
4.     Gravity  gradient  with  three  magnetic  torquers 

a.  Advantages 

Three  othogonally  oriented  magnetic  torquers  have  the  advantage  of  gen- 
erating restoring  torque  in  almost  any  direction  based  on  vector  addition.  In  this  man- 
ner, the  corrective  torque  may  be  produced  to  counter  any  disturbance  from  any 
direction  as  long  as  it  is  small  enough  not  to  tumble  the  vehicle.  Another  advantage  for 
magnetic  torquing  is  that  the  torquers  only  draw  power  while  providing  restoring  torque. 
This  amounts  to  a  considerable  power  savings  when  compared  to  the  reaction  wheel. 
The  power  drawn  by  the  magnetic  torquers  during  operations  is  based  on  the  amperes 
used.  This  allows  for  flexibility  of  design  since  doubling  the  number  of  turns  and  halving 
the  current  drawn  still  provides  the  same  magnitude  of  the  magnetic  moment. 

b.  Disadvantages 


•5' 


Magnetic  torquers  are  altitude  dependent.  The  magnitude  of  the  earth's 
magnetic  field  is  a  function  of-rr.  At  high  altitudes,  magnetic  torquing  is  not  effective. 
In  addition,  magnetic  torquing  is  not  always  available  in  all  directions.  Recalling  that 
magnetic  torque  is  Tmas  =  M  x  B,  magnetic  torque  may  not  be  generated  in  the  direction 
of  the  B  field  due  to  vector  cross  product  relationships.  Due  to  the  continuously  chang- 
ing B  field  with  respect  to  the  satellite,  this  is  a  minor  disadvantage.  The  desired  restoring 
torque  can  soon  be  generated  as  soon  as  the  B  field  direction  changes. 

5.  Thrusters 

a.  Advantages 

At  high  altitudes,  thrusters  and  momentum  exchange  devices  are  the  only 
control  methods  available  for  accurate  stabilization.  In  addition,  thrusters  are  the  most 
effective  method  for  stabilizing  a  vehicle  with  a  high  tumbling  rate. 

b.  Disadvantages 

Thrusters  are  necessary  for  satellites;  however,  any  fuel  used  for  attitude 
control  or  momentum  removal  is  less  fuel  available  for  orbit  transfer  and  maintainance. 
For  a  small,  lightweight  satellite  with  limited  fuel  storage  capacity,  fuel  resources  must 
be  economized. 

6.  Overview 

In  the  final  analysis,  the  attitude  control  system  employed  will  be  based  on 
specific  mission  requirements.  Only  when  the  ORION  mission  is  known  can  proper 
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weights  be  given  to  each  of  the  attitude  control  system  design  factors,  most  noticeably, 
power  consumption,  pointing  accuracy,  altitude,  and  desired  lifetime. 
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APPENDIX      GRADIENT  STABILIZATION  PROGRAM 

TITLE   3  AXIS   STABILIZATION  WITH  GRAVITY  GRADIENT  EFFECT 

*  DEFINITION  OF  TERMS 
* 

*  ALT=ALTITUDE  OF  SATELLITE   (METERS) 

*  R1=RADIUS   OF  EARTH   (METERS) 

*  R=RADIUS  OF  EARTH  +  SATELLITE  ALT  (METERS) 

*  RO=DENSITY  OF  AIR  AT  THE   ORBITAL  ALTITUDE 

*  WO=ANGULAR  RATE  OF  SATELLITE  WITH  RESPECT  TO  EARTH 

*  D=DRAG  FORCE=l. 7E+14*(RO/R) 

*  PSI,THETA,PHI    ARE   ANGLES   ABOUT  THE   X,Y,Z,    AXES    (RADIANS) 

*  PSIC,THETAC,PHIC   ARE   ANGLES   ABOUT  THE   X,Y,Z  AXES    (DEGREES) 

*  PSID,TD,PHID  ARE  ANGULAR  VELOCITIES  ABOUT  THE  X,Y,Z  AXES   (RADS/SEC) 

*  (B1),B(2),B(3)   ARE  THE  ANGULAR  ACCELERATION  ABOUT  THE 

*  X,Y,Z  AXES   (RADS/SEC**2) 

*  IX,IY,IZ   ARE   THE   MOMENTS   OF    INERTIA  ABOUT  THE   X,Y,Z   AXES    (KG-M**2) 

*  TX,TY,TZ   ARE   THE   SUMMATION  OF  THE   DISTURBANCE   TORQUES   AND 

*  RESTORING  TORQUES  EXCLUDING  GRAVITY  GRADIENT  EFFECTS 

*  TMX,TMY,TMZ   ARE   MAGNETIC  TORQUES   ABOUT  THE   X,Y,Z   AXES 

*  TRWX,TRWY,TRWZ   ARE   REACTION  WHEEL  TORQUES   ABOUT  THE   X,Y,Z   AXES 

*  TSX,TSY,TSZ   ARE   SOLAR   PRESSURE   TORQUES   ABOUT  THE   X,Y,Z   AXES 

*  TAX,TAY,TAZ  ARE   AERODYNAMIC   TORQUES   ABOUT  THE   X,Y,Z   AXES 

*  LX,LY,LZ    IS   THE   DISTANCE   BETWEEN  THE   CENTER  OF   PRESSURE   AND 

*  CENTER   OF  MASS    (M) 

*  MX,MY,MZ=SPACECRAFT'S   GENERATED   MAGNETIC   MOMENT  WITH  RESPECT  TO 

THE   X,Y,Z   AXES    (WB-M)      M=NIA  WHERE   M=MAGNETIC   MOMENT, 
N=#  OF  TURNS,    I=CURRENT   (AMPS),    AND  A=AREA   (M*-2).    FOR 

*  THE   MAXIMUM  VALUE   OF   M=9.0427,    N=400 , I=. 5A,RADIUS=. 12M. 

*  MG=STRENGTH  OF  EARTH'S   MAGNETIC   DIPOLE=8. OE15   WB-M  AT  R=0 

*  BX,BY,BZ=EARTH'S  MAG  FIELD  WITH  RESPECT  TO  X,Y,Z(WB/M**2) 

*  LAMS=RIGHT  ASCENSION  FROM  GREENWICH  @  TIME=0,    IN  DEGREES (LONGITUDE) 


L . ' .  .t .  jM  ju  JL  - ' -  JL  J  -  - '  -  JL  - ' -  JL  JL  JL  - '  -  - ' - -'-  - •  -  -'  - 1 


ARRAY  B(3) 

CONST  ALT=1000E+03,RO=7.0E-10,Rl=6378.E+03,.  .  . 

IX=95.  ,    IY=98.  ,    IZ=03.00,    LAMS0=0,... 

MX=9.0478,    MY=9.0427,    MZ=9.0427,    MG=8.  OE+15  , .  .  . 

KX=2.  ,    KY=2.  ,    KZ=1.  ,    TAU1=1.  ,   TAU2=1.  ,   TAU3=1.  ,.  .  . 

TSX=0.  ,   TSY=0.  ,    TSZ=0.  ,    TAX=0.  ,   TAY=0.  ,   TAZ=0. 

Jf  JLy-JLJLJLJ/.-L^.JLJLJLJ/~JLJLJLJLj*.^juJL 

*        PHI,THETA,AND  PSI   ARE   INIT  5   DEGREES  =  .08725  RADIANS  OF  ERROR* 

JL  JL  JL  J-  J-  -'-  .*-  -'-  -.'.  J/-  JL  JL  JL  JL  -'-  «L  JL  JL  JL  JL  JL  J-  JL  JL  J.  JL  JL  JL  JL  JL  JL  JL  J/-  JL  JL  JL  JL  J.  JL  «J-  JL  JL  JL  J>  JL  JL  J-  -'  -  JL  JL  J-  J-  - '  -  JL  JL  J>^  -'-  -'-  JL  JL  JL  -'-  J-  J-  JL  J,  JL  JL  JL 

INCON   ITHETA=.  08725, IPHI=. 08725, IPSI=.  08725 , ILAMS=0,PHIC=5.  ,.  .  . 
THETAC=5.  ,PSIC=5. 


DERIVATIVE 

R=R1+ALT 

D=l. 7E+14*(R0/R) 
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W0=SQRT(3.  987E+14/(R**3)) 

NOSORT 

LAMS=LAMSO+WO*TIME*(5  7.  3) 
10    IF  (LAMS. GE.  360)  THEN 

LAMS=LAMS-360. 

GO  TO  10 

END  IF 
* 
* 

*  THIS  SECTION  ADDS  THE  REACTION  WHEEL  EQUATIONS      * 

*  TRWX=-KX*(TAU1*PHID  +  PHI) 

*  TRWY=-KY*(TAU2*TD  +  THETA) 

*  TRWZ=-KZ*(TAU3*PSID  +  PS I) 

***********  jr^D  REACTION  WHEEL************************** 

******************************************************* 

*  THIS  SECTION  ADDS  THE  MAGNETIC  TORQUE  EQUATIONS 

*  AND  SOLVES  AN  EXAMPLE  FOR  THE  EARTH'S  MAGNETIC      * 

*  FIELD  (EQUATORIAL  ORBIT  ASSUMED)  * 

******************************************************* 

BX=(MG/(R**3))*SIN(.  2042)*SIN(((LAMS*PI)/180)-5. 0789) 

BY=-(MG/(R**3))*COS(.2042) 

BZ=-(MG/(R**3))*SIN(.2042)*COS(((LAMS*PI)/180)-5.0789) 


TMAGX=(MZ*BY)-(MY*BZ) 
TMAGY=( MX*BZ ) - ( MZ*BX) 
TMAGZ=(MY*BX)-(MX*BY) 

****************************** 

*  CONDITIONS  FOR  TORQUE  GENERATION  MET? 

****************************** 

IF(PHIC.GT.  .5)  THEN 

IFCPHID.  GT.  0)  THEN 

TMX=-(ABS(TMAGX)) 

ELSE 

TMX=0 

ENDIF 
ENDIF 
IF(PHIC.LT.  -.5)  THEN 

IF(PHID.LT.  0)   THEN 

TMX=ABS(TMAGX) 

ELSE 

TMX=0 

ENDIF 
ENDIF 
IFCTHETAC.GT.  .5)  THEN 

IFCTD.GT.  0)  THEN 

TMY=-(ABS(TMAGY)) 

ELSE 

TMY=0 

ENDIF 
ENDIF 
IF(THETAC.LT.  -.5)  THEN 
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IF(TD.  LT.  0)  THEN 

TMY=ABS(TMAGY) 

ELSE 

TMY=0 

ENDIF 
ENDIF 
IF(PSIC.GT.  .5)  THEN 

IF(PSID.GT.  0)  THEN 

TMZ=-(ABS(TMAGZ)) 

ELSE 

TMZ=0 

ENDIF 
ENDIF 
IF(PSIC.  LT.  -.5)  THEN 

IF(PSID. LT. 0)  THEN 

TMZ=ABS(TMAGZ) 

ELSE 

TMZ=0 

ENDIF 
ENDIF 
****Vr**ycyovvoVVnWcrt*'>v**-:VEND  MAGNETIC  TORQING  SECTION******************** 

TX=™X+TRWX+TSX+TAX 
TY=TMY+TRWY+TSY+TAY 
TZ=TMZ+TRWZ+TSZ+TAZ 

*******for  GRAVITY  GRADIENT  EFFECT  ONLY  SET  TX=TY=TZ=0************** 

B(i)=-4'V(W0-'^2)'K(IY-IX)/IZ)*PHI+W0*((IX+IZ-IY)/IX)*PSID+(TX/IX) 
B(2)  =  O*(W0**2)*((IX-IZ)/IY)*THETA+TY/IY 

B(3)=-(W0**2)^(IY-IX)/IZ)*PSI-W0--K(IX+IZ-IY)/IZ)*PHID+(TZ/IZ) 
* 

PHID=INTGRL(0. ,B(1)) 
TD=INTGRL(0. ,B(2)) 
PSID=INTGRL(0. ,B(3)) 

CPHI=INTGRL(0.  ,PHID) 
CTHETA=INTGRL(0.  ,TD) 
CPSI=INTGRL(0.  ,PSID) 
* 

PHI=IPHI+CPHI 

THETA=ITHETA+CTHETA 

PSI=IPSI+CPSI 

*PHI,THETA,PSI  ARE  =  TO  INITIAL  CONDITIONS  AS  ANNOTATED  WITH  AN  "I" 
^■-PRECEDING  VARIABLE+CALCULATED  VALUES  AS  ANNOTATED  WITH  THE  "C"  IN  FRONT 

Vc-sV  }V  VcVc  ■&  VrVr  Ve  Vr->V  Vc  -sY  >V  Vr  it  -V  it  it  VcVr  it  it  it  it  it  it  ititit  ititit  it  Vc  Vc  it  ititit  VcsYVr  it  ititit  it  A  VcVcVc  Vr  Vr-Y  }Y  ititititititititit  ititit  itit  it 

PHIC=PHI*5  7.3 
THETAC=THETA*57.3 
PSIC=PSI*57. 3 

*  USE  FOLLOWING  IF  GRAVITY  GRADIENT  EFFECT  ONLY  * 
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* 


00,DELT=200. 

,  THETAC,  PSIC, 

,  THETAC,  PSIC, 


PHIDC, 
PHIDC, 


TDC, 
TDC, 


PSIDC 
PSIDC 


****************************************************************** 

*ONTRL  FINTIM=200 
*RINT  200.  ,  PHIC 
*AVE  200.  ,  PHIC 
*RAPH  (G1,DE=TEK6 
*ABEL  (G1,DE=TEK6 
*RAPH  (G2,DE=TEK6 
*ABEL  (G2,DE=TEK6 
*RAPH  (G3,DE=TEK6 
*ABEL  (G3,DE=TEK6 
*RAPH  (G4,DE=TEK6 
*ABEL  (G4,DE=TEK6 
*RAPH  (G5,DE=TEK6 
*ABEL  (G5,DE=TEK6 
*RAPH  (G6,DE=TEK6 
*ABEL  (G6,DE=TEK6 

* 

*********************************************************** 

*    USE  THE  FOLLOWING  IF  REACTION  WHEELS  ARE  USED        * 

*********************************************************** 


18) 
18) 
18) 
18) 
18) 
18) 
18) 
18) 
18) 
18) 
18) 
18) 


TIME(UN=SEC) ,PHIC(UN=DEGREES) 

PHI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3.  00 

TIME(UN=SEC) , THETAC (UN=DEGREES) 

THETA(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3.  00 

TIME(UN=SEC) ,PSIC(UN=DEGREES) 

PSI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 

TIME( UN=SEC ) , PHIDC( UN= ' DEG/SEC  * ) 

PHID  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 

TIME(UN=SEC),TDC(UN=' DEG/SEC') 

TDC  VS  TIME  (GG^IOOOK)  IX=95  IY=98  IZ=3. 00 

TIME(UN=SEC),PSIDC(UN='DEG/SEC') 

PSIDC  VS  TIME  (GG^IOOOK)  IX=95  IY=98  IZ=3. 00 


GRAVITY  GRADIENT  ONLY  GRAPHS  ******************* 


FINTIM=400 
100.  ,  PHIC 
5. ,  PHIC, 
(G1,DE=TEK6 
(G1,DE=TEK6 
(G2,DE=TEK6 
(G2,DE=TEK6 
(G3,DE=TEK6 
(G3,DE=TEK6 
(G4,DE=TEK6 
( G4 , DE=TEK6 
(G5,DE=TEK6 
(G5,DE=TEK6 
(G6,DE=TEK6 
(G6,DE=TEK6 
(G7,DE=TEK6 
(G7,DE=TEK6 
(G8,DE=TEK6 
(G8,DE=TEK6 
(G9,DE=TEK6 
(G9,DE=TEK6 


,DELT=5. 

,  THETAC,  PSIC,  TX,  TY,  TZ 

THETAC,  PSIC,  TX,  TY,  TZ,  PHID,  TD,  PS ID 


*ONTRL 

*RINT 

*AVE 

*RAPH 

*ABEL 

*RAPH 

''"ABEL 

*RAPH 

*ABEL 

*RAPH 

*ABEL 

*RAPH 

*ABEL 

*RAPH 

*ABEL 

*RAPH 

*ABEL 

*RAPH 

*ABEL 

*RAPH 

*ABEL 
****** 

*  * 

******************************************************************** 

FOLLOWING  GRAPHS  IF  MAGNETIC  TORQUING  IS  UTILIZED  * 


TIME(UN=SEC) ,PHIC(UN=DEGREES) 
PHI(DEG)  VS  TIME  (GG@1000K)  IX 
TIME(UN=SEC) ,THETAC(UN=DEGREES 
THETA(DEG)  VS  TIME  (GG<aiOOOK) 
TIME(UN=SEC) ,PSIC(UN=DEGREES) 
PSI(DEG)  VS  TIME  (GG@1000K)  IX 
TIME(UN=SEC) ,TX(UN=NM) 
TXRW  VS  TIME  (GG@1000K)  IX=95 
TIME(UN=SEC) ,TY(UN=NM) 
TYRW  VS  TIME  (GG^IOOOK)  IX=95 
TIME(UN=SEC) ,TZ(UN=NM) 
TZRW  VS  TIME  (GG{aiOOOK)  IX=95 
TIME( UN=SEC ) , PHID( UN= ' RAD/SEC  * 
PHID  VS  TIME  (GG(aiOOOK)  IX=95 
TIME(UN=SEC) ,TD(UN='RAD/SEC' ) 
TD  VS  TIME  (GG<aiOOOK)  IX=95  IY 
TIME( UN=SEC ) , PSID( UN= ' RAD/SEC ' 
PSID  VS  TIME  (GG(aiOOOK)  IX=95 


,T.^T-,U.'«'- 


,v*********END  REACTION  WHEEL  GRAPHS***** 


=95  IY 

) 

IX=95 

=95  IY 

IY=98 

IY=98 

IY=98 

) 
IY=98 

=98  IZ 

) 

IY=98 

****** 


=98  I 

IY=98 

=98  I 

IZ=3. 

IZ=3. 

IZ=3. 

IZ=3. 

=3.00 

IZ=3. 

sV**** 


Z=3.  00 

IZ=3.  00 
Z=3.  00 
00 
00 
00 
00 


00 

****** 


*   USE 
******* 

*ONTRL 
*RINT 
*AVE 
*RAPH  ( 
*ABEL  ( 
*RAPH  ( 
*ABEL  ( 
*RAPH  ( 
*ABEL  ( 


*************Vr**Vr*************Vf*'5V*****'>V********,5'-*'>V7V********** 

FINTIM=4000.  ,DELT=5. 

100.,  PHIC,  THETAC,  PSIC,  WO ,B( 1) ,B(2) ,B(3) ,LAMS 

5. ,  PHIC,  THETAC,  PSIC,  LAMS,  TX, TY,TZ,BX, BY, BZ, PHID, TD, PSID 

G1,DE=TEK618)  TIME(UN=SEC) ,PHIC(UN=DEGREES) 

PHI(DEG)  VS  TIME  (GG@1000K)  IX=95  IY=98  IZ=3. 00 

TIME(UN=SEC) , THETAC (UN=DEGREES) 

THETA(DEG)   VS  TIME   (GG<aiOOOK)    IX=95    IY=98   IZ=3.  00 

TIME(UN=SEC) ,PSIC(UN=DEGREES) 

PSI(DEG)    VS  TIME   (GG<aiOOOK)    IX=95    IY=98   IZ=3.  00 


G1,DE=TEK618) 
G2,DE=TEK618) 
G2,DE=TEK618) 
G3,DE=TEK618) 
G3,DE=TEK618) 
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*RAPH  (G4,DE=TEK618) 
*ABEL  (G4,DE=TEK618) 
*RAPH  (G5,DE=TEK618) 
*ABEL  (G5,DE=TEK618) 
*RAPH  (G6,DE=TEK618) 
*ABEL  (G6,DE=TEK618) 
*RAPH  (G7,DE=TEK618) 
*ABEL  (G7,DE=TEK618) 
*RAPH  (G8,DE=TEK618) 
*ABEL  (G8,DE=TEK618) 
*RAPH  (G9,DE=TEK618) 
*ABEL  (G9,DE=TEK618) 
*RAPH  (G10,DE=TEK618 
*ABEL  (G10,DE=TEK618 
*RAPH  (G11,DE=TEK618 
*ABEL  (G11,DE=TEK618 
*RAPH  (G12,DE=TEK618 
*ABEL  (G12,DE=TEK618 
*RAPH  (G13,DE=TEK618 
*ABEL  (G13,DE=TEK618 
END 
STOP 


TIME(UN=SEC) 

LAMS(DEG)  VS 

TIME(UN=SEC) 

TX  VS  TIME  ( 

TIME(UN=SEC) 

TY  VS  TIME  ( 

TIME(UN=SEC) 

TZ  VS  TIME  ( 

TIME(UN=SEC) 

BX  VS  TIME  ( 

TIME(UN=SEC) 

BY  VS  TIME  ( 

)  TIME(UN=SEC 

)  BZ  VS  TIME 

)  TIME(UN=SEC 

)  PHID(RADSEC 

)  TIME(UN=SEC 

)  TD(RADSEC) 

)  TIME(UN=SEC 

)  PSID(RADSEC 


,LAMS(UN=DEGREES) 

TIME  (GG<ai000K)  IX=95  IY=98  IZ=3.  00 
,TX(UN=N-M) 
GG@1000K) 
,TY(UN=N-M) 
GG{ai000K) 
,TZ(UN=N-M) 
GG(ai000K) 
,BX(UN=W-M) 
GG@1000K) 
,BY(UN=W-M) 
GG@1000K) 
),BZ(UN=W-M) 
(GG@1000K) 

),PHID(UN='RAD/SEC') 

)   VS  TIME   (GG@1000K)    IX=95    IY=98   IZ=3. 00 
),TD(UN='RAD/SEC') 

VS  TIME   (GG@1000K)    IX=95   IY=98   IZ=3. 00 
),PSID(UN='RAD/SEC') 
)   VS  TIME   (GG@1000K)    IX=95    IY=98   IZ=3.  00 
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